Nonlinear Dynamics of Ultra-Cold Gas: Collapse of Bose Gas With Attractive 
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I. INTRODUCTION 

In recent years it has become possible to cool and to 
trap neutral atoms using the resonant laser radiation 
(see, for example; 1 and references therein). The main 
idea of cooling and trapping is connected with the de- 
pendence of radiation pressure taking into account the 
resonant absorption and photon radiation by atoms (see 
for example^). The temperature of the cooling of neu- 
tral atoms (use H and hydrogen-type atoms Na, Cs) has 
been reached up to the fraction of mK for relatively large 
samples of atoms and for minutes in a time that provided 
by engaging of atoms by magnetic trap with zero field at 
the centre of the system. The typical density n was about 



10 



li 



10 cm . For such densities the quantum effects 



become essential close to the Bose condensation temper- 
ature 



3.31S 2 n 2 / 3 
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3/2, 



and below. For example, this temperature for H has the 
order of 1 mK for such density that is more or close 
to the achieved cooling temperature. For heavy atoms 
the reached temperature until yet is far from the corre- 
sponding Tb but there exist some ideas how to get lower 
temperature^. Therefore one can hope that in the nearest 
years it will be get at the experiment the unique quantum 
gas. 

The purpose of this paper is to treat the behaviour of 
such gas for the temperature less than the condensation 
temperature when almost all particles are in the conden- 
sate state (for simplicity we restrict ourselves by T = 0) 
when the radiation is absent. It should be noted that 
resonant electromagnetic wave absorption and forthcom- 
ing reemission lead to the repulsing between atoms of 
the Coulomb type 2 . Without radiation the interaction 
between atoms in gas phase for long distance more than 
Bohr radius a is defined by the Van der Waals attraction 
forces, F — — 7/r 6 . Then in the gas approximation fol- 
lowing to^ for oscillations of the condensate one can get 



the standard Hamiltonian for the nonlinear Schrodinger 
equation (NLSE) 



H = 



2m 



|vvi s 



G 



M 4 + £M ( 



dr. 



(1) 



Here ip is the wave function of the condensate, so that 
n = \ip\ 2 is the gas density. This constant G is connected 
with the scattering length a — (m/2irh)G. The char- 
acter of the interaction depends significantly on the a 
sign. For example, for hydrogen a = 1 (in Bohr radii), 
for Cs a = 100, bu the sign is unknown. Therefore we 
consider both cases. In the repulsion case we neglect the 
next term of expansion over n proportional to f3. The 
main attention will be paid to the attraction because in 
opposite case there exists a lot of papers about, starting 
from the Bogolyubov famous paper. 



II. INSTABILITY AND COLLAPSE 

We start our consideration from the attraction case. 
The equation corresponding to H is the well-known 
NLSE, which we write in dimensionless variables for 
G < 0: 



i^ = -l/2A^-/(M 2 )V 



(2) 



where f{\ip\ 2 ) = l^l 2 — /Jl^l 4 - We shall assume that the 
initial state of the gas is homogeneous with constant den- 
sity n so that /3|V ; | 2 <C 1. Because the theory of the NLSE 
is well known for this case^r— we present here only the 
main facts from this theory paying the attention to the 
application to this concrete physical situation. 

First of all let us be reminded that such a state will be 
unstable with growth rate 



T = k(nf - fc 2 /4) 
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(3) 



(Consider tp = ip + U + iV; U, V ~ exp[rt + ikx].) The 
threshold of this instability is defined from /' = 0. It 
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gives n = 1/(2/3), that can be considered as condensed 
phase density. If we start from gas we have the instability. 
The nonlinear stage of this instability will be the collapse, 
i.e., the formation of the singularity for finite time. 

A few things should be noted about 3-D solitons as 
one possible variant of the instability development. The 
simplest solutions of this type ip = (f>exp[iX 2 i\ represent 
the stationary point of H for fixed number of particles 
N = J \^\ 2 dr: 

- A 2 + | + |^| 2 - /3|0| 4 = 0. (4) 

These solutions are unstable in gas phase according to the 
criterion: dN/dX 2 < 0^. The qualitative arguments of 
this instability can be computed from the scaling trans- 
formations, tp — > a~ 3 / 2 ip(r/a) remaining N. As a result 
of H becomes the function of parameter a: 

H = hj a 2 - I 2 /2a 3 + /3/ 3 /3a 6 (5) 

where h = ± /|W| 2 dr, J 2 = J \ip\ A dr, I 3 = J\tp\ 6 dr. 
Maximum of this function corresponds to the unstable 
soliton as to the saddle point. Besides, this function for 
f3 = occurs to be unbounded from below. It is the 
main reason of collapse in the system. From this point of 
view the collapse can be considered as the process of some 
"particle" falling in an unbounded potential. If /? is finite 
then this process can be stopped for very large amplitude 
that corresponds to the minimum of this function, i.e., 
the soliton of the condensed phase. It should be noted 
that for (3 = there is no stabilization or saturation of 
the instability. The compression will continue until the 
formation of a singularity. It follows strictly from the 
virial theorem: 

d 2 r 

— r 2 \iP\ 2 dr <4H. (6) 

After twice integration we have that positive value 
/ r 2 \ip\ 2 dr vanishes for finite time. Because for the ini- 
tial state H is negative after the formation of clusters 
with characteristic size I oc k^ 1 the collapse becomes in- 
evitable, where fco = \J2no corresponds to the maximum 
of T from Eqn. ([3]). Other questions to be raised con- 
cern the possible regimes of collapse. The first variant 
of the collapse is the so called strong regime of com- 
pression when all particles in the cavity occur in the 
singularity^. It is easy to understand that this process is 
a semi-classical one because from the usual quantum me- 
chanics it is known that in the potentials V oc — r~ a for 
a > 2 the falling becomes, near the centre, more classical. 
But such collapse is unstable with respect to short-wave 
perturbations resulting in the formation of weak singu- 
larities. For this type of collapse the captured number of 
pa! rticles is formally equal to zero, in practice it is quite 
small. The most simple argument as to why instead of 
adiabatic compression the weak singularity appears is the 
following. Let us consider some region Q with negative 



Hamiltonian H (we call this region a cavity) . Then with 
the help of the mean value theorem the inequality 

maxlVf < |flh|/2JV n (7) 

follows. Here Nq is the number of particles in SI. Let 
us imagine that from Q we have radiation of waves (or 
emission of particles). Because these emitted particles 
carry out the positive portion of H (strictly speaking it 
is valid for the separate cavity) there is a reduction in 
H for 57, i.e., H for this region becomes more negative. 
It is evident that Nq vanishes. So the ratio \Hq\/2Nq 
tends to infinity together with max \ip\ 2 - It means that 
the radiation of waves promotes the collapse, leads to its 
acceleration. The corresponding solution to this regime 
is the self-similar one of Eqn. (2) for /3 = 0: 

^ = (to - ty ll2 - lt 9 (r/(to - t) 1 ' 2 ) (8) 

where g(£) obeys the equation 

-(i-2e)g-(&)g + Ag + 2\g\ 2 g = 0. (9) 

Hence the characteristic radius of the collapsing solution 
is proportional to (to — t) 1 / 2 . From this fact follows that 
particles that reach the critical density n = 1//3 will form 
drop with the radius tq = ft 1 / 2 /\g(l)\. The resulting 
number of particles in the drop can be estimated as 

N^ \g(l)\- 3 ^Hn/3. (10) 

Because of f3n <C 1 this number is less than one. It means 
that at the collapse time to the formation of condensed 
phase does occur, secondly, the energy gathered by the 
condensate particles during the collapse is quite large and 
so there exists the unique possibility to transfer this en- 
ergy (or a significant amount of it) to undercondensate 
particles radiated from the vicinity of the singularity. It 
is very important also that while weak collapse the at- 
traction potential will be formed for other non-collapsing 
particles which can be trapped after to- It follows, in 
particular, from asymptotic of for £ — > oo: 

g(0 -+ A/e +2tt or -> \A\ 2 r- 2 . (11) 

Because during the weak collapse the number of trapped 
particles into singularity is very small, it is natural to 
suppose that the post-collapse regime is quasi-stationary. 
In searching for its distribution we can simply find a sta- 
tionary solution. It is convenient to rewrite Eqn. ([1]) for 
density n and phase $. Then simple analysis gives the 
following solution- - — : 

n(r) « l/21n(r/r )r -2 ; (12) 
v(r) = d<S>/dr « P/27rln(r/r ). (13) 

Here P is the particle flux to the centre of this black hole 
and ro the size of the growing drop. The flux, from one 
side, forms the drop and, from the other, the counter 
flux of undercondensate particles which carry out the 
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energy of the order of the kinetic energy of the falling 
particles. It is easy to estimate the corresponding en- 
ergy E lost by the cavity during the total time of col- 
lapse: E oc l-EinldSno) -1 where E m is the initial gas 
energy. Thus the carried energy appears to be quite a 
large value. It can be assumed that a sufficiently small 
portion of particles (at least, less than half of the falling 
particles) will be emitted. These particles will interact 
weakly with condensate particles. Their phases can be 
considered random and therefore their interaction will be 
described by a four- wave kinetic equation. As known 6 , 
this interaction leads to the formation of a power-type 
Kolmogorov spectrum in the particle energy space real- 
izing the particle flux in this space, i.e., forms the particle 
condensate. Now we don't know what kind of equilibrium 
is realized in this case. 



Now let us estimate the growth rate of the drop. If we 
assume that approximately half of particles from the cav- 
ity with the size I oc k^ 1 are absorbed into the singularity, 

i.e. Ndr oc tlqP oc n , then the characteristic time of 
the drop growth has the same order as the inverse growth 
rate (3): T iaax Nd r oc I 2 P. Hence P oc n^ 2 . Substituting 
this expression in the definition of P one gets the drop 

3/2 

growth as V oc (3n Q where V is the drop volume. At the 
end of this part we present the concrete numbers for an 
atom of Cs mass. For density 10 n cm~ 3 the modulation 
length / ~ 10~ 3 cm, the collapse time 0.3 sec that less 
than sample cooling size and the particle life time in the 
trap, respectively. It means that the observations of the 
discussed above effects represent quite real, of course, if 
the corresponding temperatures can be reached. 
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